Polymer transport by laminar flows 
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Polymer transport is investigated for two paradigmatic laminar flows having open and closed stream- 
lines, respectively. For both types of flows we find transport depletion owing to the action of the 
polymers elastic degree of freedom. For flows with closed streamlines the leading mechanism for the 
observed transport reduction is the (dynamical) formation of barriers. For flows with open stream- 
lines the reduction of transport is induced by the renormalization of the bare diffusion coefficient. 
Results have been obtained by means of Lagrangian simulations. 



PACS number(s) : 47.27.-i, 47.27.Qb, 47.27.Te 

The transport of particles advected by a given velocity 
field is a subject attracting increasing attention of physi- 
cists. Such issue has obvious applications in geophysics 
(e.g. pollutants dispersion) as well as in industrial pro- 
cesses (e.g. transport of a solute in porous media) Q. In 
addition, the investigation of the dynamics of particles 
leads, in a natural way, to highly non trivial behaviors as 
the anomalous dispersion and the Lagrangian chaos §§. 

In this paper we investigate the transport of passive 
polymers in given velocity fields. The interest of this 
problem is in the spectacular effect of turbulent reduction 
by the dilute polymer addities (se so-called drag reduc- 
tion 0|). Of course, the complete treatment of the drag 
reduction is very complex and must involve the back- 
reaction of the polymer concentration on the velocity field 
ruled by the Navier-Stokes equations ||U ■ 

A first crude, but non trivial, step is to consider the 
passive limit, i.e. the transport of polymers in a pre- 
scribed velocity field. This problem, in spite of its poor 
relation with the turbulent drag reduction, is not sim- 
ple at all fl§. It can be seen as a complication of the, 
already difficult, issue of the passive particle transport 

0. 

A widely used model for the polymer dynamics, due 
to Rouse [PlJlCfl, considers macromolecules as a series of 
N beads linearly connected by harmonic springs. The 
evolution law is given by the Langevin equations: 
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where V = K/2 J2n=i(^ n ~ R-n-i) 2 , Rn is the position 
of the n-th bead, u is the incompressible velocity field, 
( is the friction coefficient, K = 3ksT/b 2 is the spring 
constant, b is the average distance between the beads, ks 
is the Boltzmann constant, T is the temperature and rj n 
represents the thermal noise mimicking the interaction of 
beads with the solvent. 

The Rouse model (|l|) can be further simplified (this is 
the so-called dumbbell model) just taking N = 2: 



= K (Rr - R ) 
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We assume •q i = \/2Dof) i: with f] i normalized, zero- 
mean, independent white noise processes, and we set, 
without loss of generality, £ = 1. For a given velocity 
field, Dq and K are thus the only relevant parameters 
controling the polymer dynamics. 

In the following we shall concentrate our attention on 
the transport properties of the bead center of mass, the 
coordinates of which are denoted by R = (Ro + Ri)/2, 
R = (X, Y) . Note that in the limit of very large K the 
evolution equation for R reduces to the usual equation 
for a fluid particle but with Dq/2 instead of Dq: 



dR m\ 
-dJ = U(R) 



(4) 



As in Eqs. (g) and (g), fj is a normalized, zero-mean, 
white noise process. 

Therefore, for the understanding of the role of the elastic 
degree of freedom, results on the diffusion process of the 
bead center of mass for finite values of K have to be 
compared with those generated by the particle diffusion 
limit given by Eq. (Q). 

The model (|l|) has been investigated, e^g., in Ref. JTl] ] 
for a layered random flow, and in Ref. |12f| in the pres- 
ence of a non-potential static random flow. Antithetic 
conclusions arose for the transport properties of poly- 
mers plugged in the two above flow classes. In Ref. Jll| , 
the layered random flow has been found to cause an en- 
hanced transport with respect to the single particle dif- 
fusion problem. The opposite result (i.e., the occurrence 
of transport reduction) has been singled out in Ref. |L2| . 
For the flows considered in Ref. |12] the basic mecha- 
nism giving rise to transport reduction is played by the 
dynamically generated barriers, which are not present in 
the layered random flow |ll[] . 
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The physical key role for the barrier to emerge is orig- 
inated from the competition between flow-originated 
stretching, which acts on the polymers, and the elastic- 
ity of the polymeric structure itself. Due to the elastic 
degree of freedom, unlike particles, polymers have the 
possibility to select particular regions of the flow. Physi- 
cally speaking, polymers prefer those regions where they 
can reduce their own elastic energy. This happens, e.g., 
in those regions where the velocity strain is large and 
negative. 

The issue related to the possible existence of preferred 
regions in more realistic flows raised in Ref. |12|. One 
of the main aim here is to show, by means of numerical 
simulations, some examples of flows for which transport 
of polymers is less effective than particles transport. The 
choice for the flows falls on two paradigmatic laminar ve- 
locity fields, the diffusive properties of which have been 
analyzed in great detail in the past for what concerns 
particle dispersion |b|. To be specific, we shall focus 
on shear flows (where the streamlines are open) and cel- 
lular flows both stationary and time-dependent (where 
the streamlines are closed). As we shall see, although 
for different reasons, such classes of flow show transport 
depletion in the presence of polymers. 

Let us start from the convective flows. We investigate 
here polymers diffusion in a simple model mimicking the 
Rayleigh-Benard convection . Two-dimensional con- 
vection with rigid boundary conditions is described here 
by the following stream function: 



ij}(x,y,t) = ipo sin(a; + B sin Lot) sin y 



(5) 



where the periodicity of the cell is 2tt, and the even 
oscillatory instability jl5| is accounted for by the term 
Bsmutt, representing the lateral oscillation of the rolls. 
Velocity is obtained from (|J) by the usual relations 
u = (—dip/dy, di/j/dx). The capability of the simple flow 
(||) to capture the essential features of the convection 
problem is discussed in Ref. 14 1. 

The second flow we have considered is the two- 
dimensional Kolmogorov shear flow defined as: 



u = (u(y) , 0) 



(6) 



with u(y) — Usin(y). 

Equations (0) and @ have been integrated with u ob- 
tained from the flows (|^) and (^) using a second-order 
Runge-Kutta scheme. In what follows, averages are ex- 
tended over different realizations and are performed by 
following ~ 10 5 couple of particles. 

In order to investigate the rate of transport, a measure 
of the eddy diffusivity, D = ([X{t) - X (0)} 2 ) / (2t) , for 
large times, t, has been made, X being the x-component 
of the bead center of mass. 

Let us start the discussion of our results from the cel- 
lular flows. The behavior of the diffusion coefficient, D, 
vs the spring constant K are reported in Fig. 1 for co = 
and three different values of the molecular diffusivity. 
Dashed lines correspond to the particle limit analytically 



obtained in Ref. @ (formulae (22) and (24) with d ■ 
and = 1) in the limit of small Dq. 
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FIG. 1. The diffusion coefficient versus the spring constant 
for the convective model (||) with u = and for different 
values of Do- Dashed lines represent the theoretical values of 
the diffusion coefficient obtained in Ref. |Q for the particle 
diffusion problem. 



In Fig. 2, the behavior of the average spring elonga- 
tion, I = (|Ri — R0I 2 ) 1 / 2 is shown as a function of K, 
again for different values of Do- Dashed lines correspond 
to the limit of small I, where it is easily checked that the 
spring elongations behaves as (2Dq/ K) 1 / 2 . 




FIG. 2. The average spring elongation, I, defined as 
I = (|Ri — Ro| 2 ) 1/,2 ; versus the spring constant, for differ- 
ent values of Do. Dashed lines are relative to the limit of 
small I, when I is found to behaves as (2D0/K) 1 ' 2 . 
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From these figures, a reduction of transport always oc- 
curs for spring elongations smaller than the cell size (i.e. 
I < 7r). This is the physically interesting case. The case 
I ~ 7r corresponds to a minimum in the transport. For 
small values of Do (smaller than the smallest one shown 
in Fig. 1), there exist a narrow band of values of K cor- 
responding to I ~ 7r where particles remain trapped in 
adjacent cells. The diffusion process stops in this case. 
Note that this is a pure effect played by the elastic degree 
of freedom: for the single-particle diffusion in incompress- 
ible flows, it is possible to show that the (eddy) diffusion 
coefficient is always larger than the (bare) coefficient Do 

Transport depletion is here caused by the aforemen- 
tioned dynamically generated barriers. This can be easily 
checked in Fig. 4 where the probability density function 
(pdf ) for the bead center of mass to be in the subinterval 
[0, 7r] x [tt/4, 3/Att] of the elementary cell is shown for dif- 
ferent values of the spring constant K (integration along 
y has been performed). From this figure it appears as 
polymers prefer those regions close to the center of the 
cell where they can reduce their own elastic energy. In 
such regions, the contribution of the molecular diffusivity 
to escape from a cell to another is unimportant. On the 
contrary the regions close to the cell boundaries, where 
the role of molecular diffusivity in the (eddy) diffusion 
process is fundamental, are refused by polymers. The re- 
sult is the slowing down of the diffusing polymers. Note 
that the effect of barriers reduces as the spring constant 
increases. For K = 0.75, that corresponds to a spring 
elongation of I — 0.40, the probability density is almost 
uniform within the cell. 
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FIG. 3. The diffusion coefficient versus the spring constant 
for the 
Kolmogorov shear flow (^) for Do = 0.1. The dashed line 
represents the theoretical value, D = (Do/2) + U 2 /[2(D /2)], 
for the particle diffusion limit (see, e.g., Ref. |ff8f ). 



Behaviors similar to those shown in Figs. (1), (2) and 
(3) have been found also for lo ^= 0. As in Ref. \V7^ , we 
have restricted our attention only on the case B = 0.5. 

As one can see in Fig. 3, also in the shear case the 
diffusion coefficient, D, for finite values of K is smaller 
than that in the limit case of infinite K (see Eq. 



g 0.3 




FIG. 4. The probability density function (pdf) for the bead 
center of mass to be in the interval [0, tt\ x [tt/4, 3/47r] of the 
elementary cell. Integration along y is performed. The molec- 
ular diffusivity is Do = 0.05. For K = 0.25 we have I ~ 1.37; 
for K = 0.35, I ~ 0.75; for K = 0.75, I ~ 0.40. 



Unlike what happens in the considered cellular flows, 
this effect is not due to the presence of barriers and can 
be understood by means of simple arguments. To show 
that, introducing Sy = Y\ — Yq we immediately have from 
Eqs. (|), (| and (§): 



^ = -2KSy + y/W Q fh 



(7) 



dX 
~dt 



(u(Y - Sy/2) 



u(Y) 



u(Y + Sy/2)) 
d 2 u (Sy) 2 



dY 2 




Note that in our case = — u and, moreover, Sy per- 
forms an Ornestein-Ulenbeck process, i.e. a Gaussian 
process with (Sy) = 0, ((Sy) 2 ) = Dq/K and a character- 
istic relaxation time t — 1/(2K). Therefore, from Eq. (j|) 
one sees that the center of mass feels a renormalized ve- 
locity field u(Y) — > u(Y)(l — ((Sy) 2 )/8) and a renormal- 
ized fluctuation part, i.e. with a larger D . Now, because 
of the Taylor formula, D = D /2 + U 2 /[2(D /2)], (note 
that D — ► D / 2 in the particle limit defined by Eq. (^) 
it is easy to see that D has negative correction for finite 
values of K. 

For a generic time independent shear flow we expect the 
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same qualitative behavior. This fact can be grasped not- 
ing that for u(y) = J2 k u k exp(iky) + c.c, u k being the 
Fourier transform of u{y) and c.c. stands for the complex 
conjugate, the center of mass experiences a normalized 
velocity field £t_Ufc(l- k 2 ((Sy) 2 /8)). Exploiting the Zel- 
dovich formula f||, D ~ £ fc \u k \ 2 / (D Q k 2 ) , in all the re- 
alistic cases having rapidly vanishing values of |ufe| 2 in 
the limit k — > oo, the above argument reported for the 
Kolmogorov flow still holds. 

It is worth observing that, for time dependent shear 
flows, the increasing fluctuating part could give rise to 
transport enhancement via the interference mechanism 
identified in Ref. p0| . Indeed, in the presence of anti- 
correlated regions of the velocity field (i.e., where the ve- 
locity autocorrelation function is negative) , the enhanced 
molecular diffusivity can be advantageous to escape from 
the anticorrelated regions which slow down the diffus- 
ing particle. An enhancement in the diffusion coefficient 
might occur in this case. 

In conclusion, two different mechanisms leading to 
transport reduction have been identified for two paradig- 
matic flows with closed and open streamlines, respec- 
tively. For stationary flows with closed streamlines 
transport reduction is due to the emergence of dynami- 
cally generated barriers. For stationary flows with open 
streamlines the mechanism giving rise to the observed re- 
duction of transport is triggered by a renormalized (en- 
hanced) molecular diffusivity. In virtue of such enhanced 
molecular diffusivity, the bead center of mass forget its 
past evolution faster than in the case of particles, a fact 
that reduces the eddy diffusivity. 

It is an open question whether, in the presence of time- 
dependent shear flows, the same renormalization of the 
bare molecular diffusivity might give rise to transport en- 
hancement via the interference mechanism identified in 
Ref. @. 
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